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Abstract 

For small scattering lengths, cold bosonic atoms form a condensate the den- 
sity profile of which is smooth. With increasing scattering length, the density 
gradually acquires more and more oscillations. Finally, the number of oscilla- 
tions equals the number of bosons and the system becomes fermionized. On 
this pathway from condensation to fermionization intriguing phenomena oc- 
cur, depending on the shape of the trap. These include macroscopic fragmen- 
tation and coexistence of condensed and fermionized parts that are separated 
in space. 
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The experimental observations of Bose- Einstein condensates (BECs) [1,2] have promoted 
a contemporary search for the physics governing trapped cold bosonic systems [3] . For weakly- 
interacting bosons, Gross-Pitaevskii theory [4] provides an excellent description of the ground 
state. For the other extreme case, namely, for impenetrable bosons in one dimension (ID) - 
so-called Tonks-Girardeau gas - the mapping between bosons and free (spinless) fermions is 
exact [7] and provides the ground state. Clearly, the Gross-Pitaevskii theory must fail in ID 
somewhere in between these extreme cases, also see [5,6]. The crossover from a BEG to the 
Tonks-Gerardeau regime has been a subject of considerable interest and its study has mainly 
focused on harmonic traps [6,8-13]. With the recent realization of a Tonks-Girardeau gas 
[14], it is anticipated that theoretical predictions would be captured in the lab. 

What happens on the pathway from condensation to fermionization and how, at all, to 
follow this pathway for a general trap potential have yet to be addressed. In fact, even for 
harmonic traps oscillations of the density due to "fermionization" of the bosonic system have 
not been recovered for finite scattering lengths. This Letter addresses the crossover from 
condensation to fermionization in a general trap potential. We will show that (a) the density 
gradually acquires oscillations as the scattering length increases, the number of which equals 
the number of bosons when the system becomes fermionized (enters the Tonks-Girardeau 
regime) ; (b) the pathway from condensation to fermionization depends on the shape of the 
potential. In a general trap, intriguing phenomena not encountered in the homogeneous 
or harmonic system can occur. For instance, in asymmetric double-well potentials a cold 
bosonic system can be condensed (coherent) in one well and fermionized in the second one. 

We would like to follow the pathway from weakly- to strongly-interacting bosons in real 
space and obtain a wavefunction picture of the ground state of cold atoms in the trap. This 
will allow us to monitor directly the fine changes in the spatial density of interacting bosons 
as the inter-particle interacting increases. Gross-Pitaevskii theory, which is an excellent 
mean-field for weakly-interacting boson systems, fails to describe the system as the interac- 
tion increases [5,6]. There is, however, a more general mean-field for bosons that employs a 
single configurational wavefunction (also termed permanent) [15]. Generally, we may put rii 
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bosons to reside in one orbital, 0i(r), n2 bosons to reside in a second orbital, 02 (r), and so 
on, distributing the atoms among 1 < Uorb < N orthogonal orbitals. We are going to look 
for the ground state of the iV-boson system in the subspace of all possible configurations: 
* = <S0i(ri) • • • 0i(r„j02(rni+i) • • • hir "1+712) ■ ■ ■ 0no^i,(rAr), where S is the symmetrization 
operator. Gross-Pitaevskii theory - where all bosons reside in a single orbital - is now only 
one possible configuration in which the system might be found. 

Let the many-body Hamiltonian describing N bosons in a general trap potential be 
H — J2f=i T{ri) + V{ri) + Y.ii.j=i U (r^ — r^). Here, rj is the coordinate of the i-th particle, 
T{ri) and V{ri) stand for the kinetic energy and trap potential, respectively, and U (r^ — rj) = 
gS{ri — Tj) describes the pairwise contact interaction between the i-th and j-th atoms, 
where g is proportional to the s-wave scattering length. Taking the expectation value of the 
Hamiltonian with respect to ^ leads to the following energy functional [15]: 

EN[norb,{n^},{<P^}] = T.[^^ J ^K^) [t {r.) + V 0.(r)cir + ^^^^^^^^ J |0.(r)|^rfr| 

+ J22gn,nJ IMrmM'dr. (1) 

The A'"-boson energy functional depends on the number riorb of orthogonal orbitals 
in which the bosons reside, the occupations {n,} of these orbitals and the orbitals 0j(r) 
themselves which have to be determined self- consistently. In order to find the ground state 
of the many-bosonic system in the subspace of all possible configurations one has to 
minimize the energy functional (1) with respect to its arguments. In practice, this results in 
Uorh coupled, non-linear equations which have to be solved self- consistently. The occupations 
{ui} are determined variationally to minimize the energy. 

The multi-orbital energy functional (1) has been successfully employed and led to the pre- 
diction of (macroscopic) fragmentation of the ground state of repulsive condensates [16] and 
of additional quantum phases of cold bosonic atoms in optical lattices [17]. In macroscopic 
fragmentation, a large number of atoms reside in a small number of orbitals. Specifically, 
three orbitals were considered in [16] within Eq. (1). Wc shall see below that, on the pathway 
from condensation to fermionization the number of orbitals needed to describe the ground 



state gradually grows. 

In what follows, we study the "evolution" of a trapped bosonic system as the inter- 
action strength increases. This will allow us to explore the pathway to fermionization in 
detail. We will see that in between condensation and fermionization the shape of the con- 
fining potential is particularly important. To demonstrate and explore the reachness of 
possibilities for a bosonic system before it fermionizes, we chose as an illustrative ID ex- 
ample the potential profile V{x) shown in Fig. 1. It is an asymmetric double- well potential 
with the left well being deeper and narrower than the right one. The problem is rescaled 
to the dimensionless coordinate x in which the ID many-body Hamiltonian is written as 
H = Ef [-d'/dx^ + V{xi)] + Ei^,- Xo5{xi - xj). 

For transparency of presentation, in the numerical example presented we follow the 
pathway from condensation to fermionization of iV = 12 bosons. We have obtained similar 
results for various numbers of bosons. We begin with a small value of the inter-particle 
interaction Aq. For convenience we introduce the quantity A — '^"^^""'"^ ■ For A = 1, the 
energy functional (1) is minimal when all atoms reside in one orbital, thus recovering Gross- 
Pitaevskii theory. The bosonic system is a condensate which resides in the left, deeper 
potential well, see Fig. 1(a). If we enforce the bosons to occupy more orbitals, the energy is 
found to increase. In other words, for A = 1 the best mean-field is Gross-Pitaevskii theory. 
It is instructive to briefly analyze the energy functional E^, Eq. (1), for a general system of 
weakly-interacting trapped bosons. For weakly-interacting bosons, the one-particle terms 
dominate E^. Hence, bosons prefer to accumulate in one orbital which is lowest in energy. 
This is exactly the Gross-Pitaevskii energy and the energy functional (1) boils down to it in 
case of weakly-interacting bosons. 

Next, when the inter-particle interaction is increased beyond A ~ 18 our system becomes 
two-fold fragmented. We remind that the fragmentation of a repulsive condensate in the 
ground state was discovered only recently [16], and is by itself a basic phenomena to be 
further explored. For A = 20, ten bosons occupy one orbital which is localized in the 
left, deeper well, whereas 2 bosons occupy now an orthogonal orbital localized in the right, 
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shallower well, see Fig. 1(b). For A = 30 the system remains two- fold fragmented with 8 
bosons in the left well (orbital) and 4 bosons in right well (orbital). For A = 40 and 50 
the two-fold fragmented ground state consists of 7 bosons in the left well and 5 bosons in 
the right well, see Fig. 1(c). So far, in the regime up to A = 50 the system is a two-fold 
fragmented condensate. Each of the two orbitals occupied by the bosons is mainly located 
in one of the two distinct wells of the potential. If we enforce the bosonic system to occupy 
more orbitals, the energy becomes higher. This means that the best mean-field here is a 
two orbital mean- field, which we denote as MF(2). In other words, Ejsr obtains its minimum 
with two self-consistent orbitals and their corresponding occupations. 

Let us briefly explain why fragmentation of the BEC occurs in the double-well potential. 
The mean-field boson- repulsion term in the Gross-Pitaevskii energy functional [Eq. (1) with 
Uorb = 1] scales with the number of particles as Ao-iV(A^ — 1)/2. When the repulsion between 
the bosons increases, the energy quickly grows. If, however, we distribute the bosons between 
two orbitals localized in the left and right wells with occupations of rii and n2 — N — rii, the 
main contribution from mean-field terms in (1) scales as Aq • [N{N — l)/2 — ani{N — rii)], 
a ~ 1. From a certain value of the interaction strength on, it becomes preferable to occupy a 
second orbital since the cost in kinetic energy is compensated by the smaller boson-repulsion 
energy. 

What happens to the bosonic system when A is further increased? With increasing inter- 
action, the distribution of the bosons between two fragments (orbitals) becomes energetically 
unfavorable, and a third orbital starts to be occupied, see Fig. 1(d) for A = 60. There are 
still 5 and 7 bosons in the left and right wells, respectively, as found in the bifragmented 
condensate for A = 40 shown in Fig. 1(c), but the bosons in the right well occupy now two 
orbitals with the occupations 6 and 1. When A is further increased, the fragment in the 
right well further breaks apart, whereas the fragment in the left well remains unchanged, 
see Fig. l(e,f). For A = 100, almost all bosons in the right well reside in different orbitals, 
whereas all bosons in the left well still reside in a single orbital. This ground state of the 
system is very surprising. It shows the coexistence of a coherent part in the left well and a 
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(almost) fully fermionized part in the right well. 

We can explain this intriguing situation that the trapped boson system is found in 
a ground state with coexisting condensed and fragmented parts as described above, see 
Fig. 1(f). For this, we remind that the two fragments of the double- well system are nearly 
non-overlapping and hence can be crudely viewed as two "almost independent" condensates. 
The effective strength 7 (see below) of the inter-particle interaction in the right well is larger 
than that in the left well, because the density in the right well is smaller. Hence, the right 
fragment almost finishes "fermionizing" while the fragment in the left well still remains 
"condensed" . 

For A = 110, each boson in the right well occupies a different, orthogonal orbital. On the 
other end, the fragment in the left well just begins to spilt, with 4 and 1 bosons residing in 
two different orbitals, see Fig. 1(g). For A = 120, one boson moves from the left to the right 
well which now houses 8 bosons in 8 different orbitals, and in the left well there is a 3-boson 
pack and a single additional boson, see Fig. 1(h). For A = 150, we find in the left well a 
2-boson pack. All remaining bosons, 2 more in the left well and 8 in the right well, reside 
in different orbitals, see Fig. l(i). Finally, from A ^ 155 on, all bosons reside in different 
orthogonal orbitals - the bosonic system has become fully fermionized, see Fig. l(j). It is 
instructive to express this critical value of A in terms of the dimensionless parameter 7 often 
employed in ID [6] which reads here 7 = Ao/2n, where n is the boson density. When this is 
done, we obtain for the present bosonic system that onset of fermionization is at 7 ~ 9, a 
value that can be obtained nowadays experimentally (see Ref. [14]). 

When more particles are trapped in the double- well trap employed here, the physical pic- 
ture obtained above remains essentially the same - we checked it numerically. The pathway 
we obtained starts, of course, from condensation, passes through fragmentation, gradual 
build-up of fermionization and eventually arrives at (full) fermionization. The regime of 
coexistence of fully condensed and fermionized parts is relatively sensitive to the number of 
trapped bosons and also to the shape of the trap. We note in this respect that a coexistence 
of ui condensed bosons and N — ui fermionized bosons can be realized by appropriately 
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designing double-well traps for any numbers N and rii < N oi atoms. 

The description of a bosonic system as gradually occupying more orbitals with increasing 
scattering length, all the way from ideal condensation {uorb = 1) to fermionization {uorb = 
N) is physically very appealing. It allows us to microscopically monitor the competition 
between one-body kinetic and potential energies and the two-body repulsion in a genera/ trap 
potential. What is important to discuss is how this pathway may be observed experimentally. 
Shown in Fig. 2 is the density profile of the bosonic system, p{x) — jj: 'ni\(f)i{x)\^ which 
can be measured. As the interaction strength increases and more orbitals becomes occupied, 
the density gradually acquires more oscillations in real space. The number of oscillations 
increases with increasing A, see Fig. 2, until it saturates at N - the number of bosons in 
the system. This, as just stated above, is when the system becomes fermionized. Another 
possibility to follow experimentally the pathway to fermionization could be to release the 
trap and register the interference patters obtained for different A's which should gradually 
loose contrast as A grows. 

The pathway from condensation to fermionization can be obtained for any bosonic system 
by minimizing the energy functional (1). Let us examine the pathway in terms of the energy 
functional E^. As the scattering length increases, the two-body, mean-field energy terms 
start to dominate En and the occupations {rii} of the bosons would like to become smaller 
in order to minimize the energy. This leads to the occupation of more and more orbitals 
in the ground state. Eventually, all rii — 1 and the last term in the braces of Eq. (1) 
becomes identically zero, thus reversing the situation described by Gross-Pitaevskii theory 
where all bosons reside in a single orbital. When this happens, i.e., when each boson resides 
in a different, orthogonal orbital, the system becomes fermionized. What happens when 
the inter-particle interaction is further increased? As this happens, the Uorb = N occupied 
orbitals become more localized in the sense that their overlaps (the orbitals seen in Fig. 1 
overlap and hence exhibit negative values as well) are gradually reduced. This prevents the 
last term in (1) from diverging. At the extreme case of infinitely-strong interaction there 
is no overlap at all between the orbitals, the last term in (1) drops out and the energy 

7 



saturates. This generic "evolution" from weakly- to strongly-interacting bosons holds true 
for any trap potential and any dimension, what makes the energy functional (1) promising 
also beyond ID. 

In conclusion, in this work the evolution of the density profile of a cold bosonic system 
with increasing scattering length is followed. With increasing scattering length the density 
profile acquires more and more oscillations, until their number eventually equals to N , the 
number of bosons. The ground state and density profile of a bosonic system strongly depend 
on the shape of the trap potential. This allows one to design systems with an intriguing 
ground state, e.g., for which one part is condensed and the other one is fermionized. The 
employment of a more general self- consistent mean-field theory has been shown to be able 
to follow - in real space - the "process" of fermionization of a general cold bosonic system. 
As to experimental realizations: Double-well potentials can be prepared by optical dipole 
traps [18] or by combining an optical lattice and strong harmonic confinement [19]. The 
scattering length can be increased then by Feshbach-resonance techniques [20] . Alternatively, 
micorfabricated traps - so-called atomchips - offer a promising venue to more general low- 
dimensional trap potentials along with controllability of the effective scattering length over 
a wide range of interaction strengths [21,22]. 
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FIG. 1. Ground-state wavefunction (orbitals (j)i{x) and occupations) of = 12 bosons in an 
asymmetric double-well potential V{x) (scaled here by 1/100) for interaction strengths: (a) A = 1; 
(b) A = 20; (c) A = 40; (d) A = 60; (e) A = 80; (f) A = 100; (g) A = 110; (h) A = 120; (i) 
A = 150; (j) A = 160. As A increases, additional orbitals are needed to describe the ground state 
which turns from a condensate to a fragmented condensate and eventually to a fermionized state. 



For convenience, the \Jni/N(j)i{x) (scaled here by 1/3 and shifted vertically and equidistantly) are 
plotted. The energy-per-particle e for each A is indicated on the plot. 
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FIG. 2. Ground-state density p{x) = ^2^°'''' of N = 12 bosons in the asymmetric 
double-well potential V{x) shown in Fig. 1 for interaction strengths in the range A = 1 to 160. For 
convenience, p{x) is scaled here by 1/2 and shifted vertically. 
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